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ABSTRACT. Population dynamics of the epiphytic bromeliad, Tillandsia brachycaulos, was followed for a 
year (October 1996 to October 1997) in a low-statured tropical deciduous forest in the Yucatan Peninsula, 
Mexico. Nine hundred and twenty individuals of this species, located on nine trees of Gymnopodium 
floribundum (Polygonaceae), were tagged. The individuals were classified into two major groups according 
to their origin: seed or vegetative propagation. Each group was subdivided with regard to height. Bimonthly 
demography censuses were carried out, where growth, survival, and fecundity were reported for each 
individual. Results indicate that most of the individuals arise from vegetative propagation (66%), highest 
mortality value was found in the seedling stage (71%), highest fecundity in the adult 3 stage (2.05 seedlings/ 
individual), whereas the highest production of offshoots is presented by the adult 3 offshoot stage (1.13 
offshoots/individual). A Lefkovitch matrix model was used to estimate the population finite rate of growth 


(A), that resulted in a value of 0.91. 
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INTRODUCTION 


In the last few decades, the accelerated loss 
of plant and animal species caused by human 
activities, has reached alarming numbers. It is 
estimated that 23% of all vascular plant species 
will have been lost by the end of the 20th cen- 
tury (Raven 1988). Some plant families such as 
the Orchidaceae, Bromeliaceae, Arecaceae, and 
Cactaceae have been threatened for some time. 
Therefore, one of the human priorities must be 
the conservation of species (Ehrlich 1988, Nad- 
karni 1992, Primack 1993). 

Schemske et al. (1994) suggested that the re- 
covery efforts for an endangered species require 
answering the following three questions: What 
is the biological status of the species? What are 
the life cycle stages that have the greatest effect 
on population growth and species persistence? 
What are the biological causes of variation in 
those life history stages that have a major de- 
mographic impact? Two of these questions can 
be answered using a demographic approach and 
incorporating the projection of population ma- 
trix models. 

Projection matrix models can be used to ex- 
plore the biological status of a population and to 
identify the life history stages that most affect 
the population finite rate of increase (A) (Cas- 
well 1989, van Groenendael et al. 1988). Tran- 
sition matrices are used to estimate à while elas- 
ticity matrices are used to measure the change 
in the value of à in response to relative changes 
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in the value of a matrix element (de Kroon et 
al. 1986, Caswell 1989). 

Some population dynamics studies that use 
projection matrix models attempt to estimate the 
regenerative capacity of the species and deter- 
mine which are the most vulnerable stages of the 
life cycle, with the objective of applying this 
information to conservation and management 
plans for several species (Peters 1991, Pinard & 
Putz 1992, Durán & Franco 1995, Martinez & 
Alvarez-Buylla 1995, Olmsted & Alvarez-Buyl- 
la 1995). 

Unfortunately, there are few studies of popu- 
lation dynamics of epiphyte species (Benzing 
1981, Larson 1992, Hietz 1997), and even fewer 
that use the projection matrix models (Tremblay 
1997, Hernandez-Apolinar 1992). Epiphytes are 
organisms that spend a major portion of their life 
cycle situated on other plants without taking nu- 
trients directly from the host plant (Benzing 
1990). They play an important role in tropical 
ecosystems where they contribute significantly 
to the richness and abundance of species (Gen- 
try & Dodson 1987). 

Habitat destruction and illegal exploitation 
have threatened many populations of epiphyte 
species (Hernández-Apolinar 1992, Larson 
1992, Rauh 1992). In consideration of this sit- 
uation, the generation of information that helps 
to understand the behavior of the populations of 
epiphytic species is of crucial importance. 

Among vascular plants, Orchidaceae, Ara- 
ceae, and Bromeliaceae contain the highest per- 
centage of epiphytic species (Kress 1986). In the 
Bromeliaceae, the neotropical genus Tillandsia 
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comprises the largest number of epiphyte spe- 
cies. Unfortunately, many species of Tillandsia 
have become fashionable ornamental plants on 
the European and Japanese markets (Rauh 
1992), which has occasioned the illegal collec- 
tion of many species. Seven species of Tilland- 
sia are listed in CITES as endangered and 14 
more species have been suggested by Rauh 
(1992) to be included. 

The purpose of this study is to analyze the 
population dynamics of Tillandsia brachycau- 
los, determine the status of populations in Dzi- 
bilchaltun National Park, and learn which are the 
most vulnerable stages of the life cycle. 


MATERIALS AND METHODS 


Fieldwork took place from October 1996 to 
October 1997 in a low-statured dry tropical for- 
est in the National Park of Dzibilchaltün in the 
state of Yucatán in Mexico (21?05'N, 89?93'W). 
According to García (1964), the study site is in 
an area with a subhumid climate, with a dry sea- 
son usually lasting from December to May. Av- 
erage annual precipitation is 879 mm, average 
annual temperature is 25.8°C (Thien et al. 1982). 
The vegetation is deciduous, with a maximum 
canopy height of 8 m. Bursera simaruba (Bur- 
seraceae), Caesalpinia gaumeri (Caesalpini- 
aceae), Piscidia piscipula (Fabaceae), and Gym- 
nopodium floribundum (Polygonaceae) are com- 
mon tree species. The epiphytic community is 
formed by species of Cactaceae, Orchidaceae, 
and Bromeliaceae, the latter being one of the 
most diverse families in the area. Within the 
Bromeliaceae, the genus Tillandsia is represent- 
ed with the most species: 7. paucifolia, T. re- 
curvata, T. schiedeana, T. elongata var. imbri- 
cata, T. fasciculata, and T. brachycaulos, the lat- 
ter being the most abundant at the study area 
(pers. obs.). 

Tillandsia brachycaulos is an epiphytic spe- 
cies, native to Southern Mexico, and Central 
America to Panama. The plant is a stemless 
monocarpic rosette with numerous subrosulate 
12-26 cm long leaves, the laminae densely and 
finely lepidote. All leaves turn an attractive red- 
dish crimson color at flowering time, and from 
the center of the rosette emerge a dozen erect 
flowers with linear, lilac-colored, tubular petals 
5—7 cm long; the stamens and pistil are exerted; 
the capsule is cylindrical, acute, 4 cm long. Af- 
ter flowering, the rosette dies and produces 1 or 
2 offshoots (pups) in the axil of basal leaves 
(Utley 1994). In Dzibilchaltun, 7. brachycaulos 
is found growing on phorophytes such as Phyl- 
lostylon brasiliense (Ulmaceae), and Acacia 
gaumeri and Pithecellobium dulce (Mimosa- 
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TABLE 1. Categories (life cycle stage) and their size 
ranges in Tillandsia brachycaulos populations in 
Dzibilchaltun National Park based on rosette 
height. 


Height range 


Category (cm) 
Seedling (s) >0.4 
Infantile (1) 0.5-2.4 
Juvenile (j) 2.5—4.0 
Adult 1 (a1) 5-7.9 
Adult 2 (a2) 8—11.4 
Adult 3 (a3) 11.5 
Infantile offshoot (10) 0.5-2.4 
Juvenile offshoot (jo) 2.5—4.0 
Adult 1 offshoot (aol) 5—7.9 
Adult 2 offshoot (ao2) 8—11.4 
Adult 3 offshoot (ao3) plants >11.5 


ceae), but is most commonly found on Gymno- 
podium floribundum. 

In order to understand the population dynam- 
ics of Tillandsia brachycaulos in the park, fe- 
cundity, survival, and growth of all accessible 
individuals of this species, located on nine trees 
of Gymnopodium floribundum (in three plots 
randomly selected to cover spatial variability), 
were followed for a period of a year. Each in- 
dividual was tagged with a number. The individ- 
uals were classified into two major groups ac- 
cording to their origin, i.e., individuals arising 
from seeds or individuals arising by vegetative 
propagation. Each group was subdivided into 
height categories (measured from the base of the 
rosette to the top of it). Each category within 
each group represents individuals with similar 
growth rates, survival, fecundity, and vegetative 
propagation. In order to compare between 
groups, analogous stages were created within 
each one, except for the seedling stage. Thus, 
plants were placed in one of the eleven stages 
(TABLE 1) and the population structure was de- 
termined. 

A census of all tagged individuals was carried 
out every two months where height, rosette ar- 
chitecture, survival, number of capsules, and 
offshoots produced by each individual were 
scored. 

With the information of the census, the fate 
of each individual after a year was determined. 
The individuals may have different fates. They 
may transit from one stage to another, stay in 
the same stage, die, or produce seeds and/or off- 
shoots. With these data, the different probabili- 
ties for the population were performed: transi- 
tion (g) was calculated as the proportion of in- 
dividuals of each stage that moved to another 
stage; stasis (s) was estimated as the proportion 
of individuals that stayed alive in the same 
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FIGURE 1. Population structure in the October 1996 


census. In each category the first column represents 
the individuals originated from seed, the second col- 
umn represents individuals originated from vegetative 
propagation. Categories as in TABLE 1. 


stage. The fecundity (F) was calculated on the 
basis of the average number of seeds/individual 
per stage, multiplied by 0.25. Based on field ob- 
servations, this value was arbitrarily assigned, 
and represents a very optimistic estimation with 
regard to the dispersal of seeds to a safe site (to 
have a better estimate of this probability, a seed 
dispersal experiment is being carried out). The 
resulting number was multiplied by 0.024. This 
value represents the percentage of germination 
obtained in a field experiment. Vegetative prop- 
agation per stage (f) was estimated as the num- 
ber of offshoots per individual. 

Using field data, a life-cycle diagram was cre- 
ated and the values of the demographic process- 
es (growth, stasis, fecundity, and propagation) 
were used to make a Lefkovitch matrix model 
to estimate the finite growth rate of the popula- 
tion (A). | 


RESULTS 


The original population structure, established 
in October 1996 using the categories shown in 
TABLE 1, is represented in FIGURE 1. Most of the 
individuals originated from vegetative propaga- 
tion (66%), and 34% of the individuals origi- 
nated from seed. Considering both groups, the 
intermediate stages (juvenile, adult 1, adult 2) 
are the most abundant. The latter three catego- 
ries, resulting from vegetative propagation 
alone, compared to the total number of individ- 
uals of the population (920) made up 51% of the 
individuals. Twenty-three percent of the individ- 
uals originated from seed. Infantile stages were 
well represented (6.6% and 8.4%) as well. The 
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FIGURE 2. Life-cycle diagram of Tillandsia bra- 
chycaulos. Only adult individuals introduce new indi- 
viduals to the population, through seed and/or vege- 
tative propagation: g (solid line) indicates the proba- 
bilities of moving from one stage to another; s (dotted 
line) represents the probabilities of staying in the same 
stage; F (long dash/dotted line), the number of seed- 
lings for each individual sampled stage; f (dashed line) 
represents the number of offshoots for individuals on 
each sampled stage. Categories as in TABLE 1. See text 
for details. 


least represented is the seedling stage (0.15%) 
followed by adult 3 stage (1.9%). 

The estimated values for mortality, fecundity, 
and propagation are displayed in TABLE 2. Based 
on the monocarpic nature of the rosettes, the 
highest value of mortality corresponds to cate- 
gories adult 2, adult 3, adult 2 offshoot, and 
adult 3 offshoot. Among the non-adult stages, 
the highest mortality is found in the seedling 
stage (71%), where most of them die when they 
fall off the phorophyte. 

The highest values of fecundity and vegeta- 
tive propagation are presented in the largest 
stages (a2, a3, ao2, ao3). Even though the pro- 
duction of offshoots normally occurs after the 
individual has flowered, 1t may also occur when 
individuals die without flowering, even in the 
infantile stage. 

The life cycle diagram of Tillandsia brachy- 
caulos is shown in FIGURE 2. It is a complex 
diagram which indicates the fate of individuals 
originating from two different forms of regen- 
eration, vegetative propagation, and reproduc- 
tion by seeds. The diagram shows the possible 
transitions that may occur with the individuals 


Values of mortality (76 M), fecundity (E as no. of seedlings/individual), and vegetative propagation 


(P) (No. offshoots/individual), as occured from October 1996 to October 1997. Categories as in TABLE 1. 


Stages S i J al a2 
M 71 40 34 69 85 
F 0 0 0 0.394 0.94] 
P 0 0.039 0.108 0.448 0.754 


a3 io jo aol ao2 ao3 
100 22 43 60 88 100 

2.02 0 0.05 0.596 1.369 1.966 

0.94 0.246 0.3 0.462 0.897 1.131 
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Projection matrix corresponding to the life cycle diagram of Tillandsia brachycaulos. Categories as _ 


in TABLE 1; N represents the number of individuals per category found in October of 1996. 


S 1 ] al a2 
S 0 0 0 0.394 0.941 
1 0.286 0.167 0 0 0 
j 0 0.33 0.189 0 0 
al 0 0.103 0.378 0.134 0 
a2 0 0 0.081 0.179 0.066 
a3 0 0 0.014 0 0.082 
10 0 0.039 0.108 0.448 0.754 
jo 0 0 0 0 0 
aol 0 0 0 0 0 
ao2 0 0 0 0 0 
ao3 0 0 0 0 0 
N 14 78 74 67 61 


of Tillandsia brachycaulos during one year. 
Along the solid arrows (g) an individual may 
advance more than one category, especially in- 
dividuals originated from vegetative propaga- 
tion, or may stay in the same category (s). At 
the same time, one may observe that an adult 
individual may incorporate new individuals ei- 
ther as seedling (F) or as offshoot (f). 

The analysis of the transition matrix (TABLE 
3) resulted in a finite growth rate (A) of 0.91. 
The stable distribution (w) (FIGURE 3), that rep- 
resents the proportional distribution of individ- 
uals of each stage when the value of à is stable, 
shows that the population of Tillandsia brachy- 
caulos will be made up mostly of seedling 
(30%) and infantile offshoot (29%) stages. The 
high number of seedlings calculated with the 
stable population structure (w) was strongly in- 
fluenced by the probability of dispersion (0.25) 
used in the calculation of the fecundities. If the 
probability of dispersion factor were to be re- 
duced, the stable population structure (w) would 
change. The proportion of the individuals of the 
rest of the stages is more or less similar for in- 
dividuals originated from vegetative propagation 
and for individuals originated by seeds. 
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FIGURE 3. Stable structure of the population (w). 


The first column represents the proportion of individ- 
uals of the population originating from seed, the sec- 
ond column represents the proportion of individuals of 
the population originating from vegetative propaga- 
tion. Categories as in TABLE 1. 


m 


a3 iO jo aol ao2 ao3 
2.05 0 0.052 0.596 1.369 1.96 
0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 
0.94 0.246 0.3 0.462 0.897 1.13 
0 0.033 0.57 0 0 0 

0 0.148 0.183 0.035 0 0 

0 0.18 0.207 0.193 0.037 0 

0 0.115 0.126 0.175 0.066 0 

8 61 87 171 213 76 


The reproductive value at stable distribution 
(v) measures the value of an individual as a seed 
for future population growth, and represents the 
amount of future reproduction, the probability of 
surviving to realize it, and the time required for 
the offspring to be produced (Caswell 1989). As 
is shown in FIGURE 4, the highest reproductive 
values for the population of Tillandsia brachy- 
caulos are in the oldest categories, in both types 
of individuals. 


DISCUSSION 


This study presents the preliminary results of 
the first study of the population dynamics of a 
species of Tillandsia using matrix models. Til- 
landsia brachycaulos, even though it is not con- 
sidered to be endangered, is a commercial spe- 
cies and usually abundant, the latter being the 
reason why it is extracted from natural popula- 
tions. 

This study shows some interesting results. To 
begin with, the original population structure dis- 
plays two patterns: one for individuals originat- 
ing from seed having a descending pattern, i.e., 
the largest numbers of individuals are in the 


0.15 - 


0.05 


Reproductive value (v) 


Ho  a1-ao1 a2-ao02 a3-ao3 


Category 


I-lo 


FIGURE 4. Reproductive value (v). The first column 
represents individuals originating from seed, the sec- 
ond column represents individuals originating from 
vegetative propagation. Categories as in TABLE 1. 
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younger categories (except seedlings). The other 
pattern, i.e., of individuals arising from vegeta- 
tive propagation, presents the majority of the in- 
dividuals in the intermediate categories. This 
suggests that the selection pressures on the in- 
dividuals of the two origins of these populations 
are very different. 

The population structure encountered in Oc- 
tober of 1996 is very different from the stable 
stage structure predicted by the analysis. This 
suggests that the population is far from being in 
equilibrium. The latter may be due to the fact 
that the rate of growth, fecundity and mortality 
Observed were different from those of earlier 
years, probably due to the arbitrary factor of dis- 
persion (0.25) used in the analysis, or that per- 
turbations occurred which changed the popula- 
tion (Mendoza 1994, Martínez & Alvarez-Buyl- 
la 1995). 

Another important aspect is that vegetative 
propagation plays an important role in maintain- 
ing the population of Tillandsia brachycaulos. 
The high proportion of individuals originating 
from vegetative propagation sustains this state- 
ment. Contrary to data reported for other species 
with vegetative propagation (Lovett 1981, Aug- 
spurger 1985, Mendoza 1994, García-Franco & 
Rico-Gray 1995), the recruitment of new indi- 
viduals originating from seed in 7. brachycaulos 
is low, but frequent. These findings agree with 
Benzing's (1981) suggestion that epiphytic 
plants require an extraordinary combination of 
vegetative propagation (to secure colonized 
sites) and the production of vagile seeds to con- 
quer new sites. 

The few seedlings found, and their high mor- 
tality, suggest that the bottleneck of Tillandsia 
brachycaulos populations is the establishment 
and survival of the seedlings, as in many other 
epiphytic species. Benzing (1978), in a germi- 
nation experiment of Tillandsia circinnata (Bro- 
meliaceae), found that although germination in 
the laboratory was high, in the field it only 
reached 3.04%, and the seedlings had low sur- 
vival in the field (Benzing 1981). Hernández- 
Apolinar (1992) reported similar results with 
germination of 100% in laboratory and only 3.5 
X 10% in the field for Laelia speciosa (Or- 
chidaceae), while Larson (1992) reports low 
seedling survivorship (15.5%) for Encyclia tam- 
pensis (Orchidaceae). Ackerman and Zimmer- 
man (1996) report the same tendency for Tol- 
umnia variegata (Orchidaceae). 

The finite population growth rate (A = 0.91) 
indicates that the population is decreasing. Giv- 
en the abundance of individuals of Tillandsia 
brachycaulos which is observed in the Park, one 
would expect that the population should be in- 
creasing or at least be at a stable stage. The rea- 
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son that the value of lambda is less than 1 could 
be due to the fact that the mortality rates were 
very high (0.67), which may be the result of the 
fact that the majority of the individuals (66%) 
were found in the reproductive stages (al, a2, 
a3, aol, ao2, and ao3). Given the monocarpic 
nature of the rosettes, they die shortly after fruit- 
ing. The high percentage of fallen individuals 
(11.3%) this year may also be responsible. 

The preliminary results of this study help us 
in answering two of the questions made by 
Schemske et al. (1994), for the implementation 
of management plans and conservation of a spe- 
cies. We can say that the population of Tilland- 
sia brachycaulos in the National Park of Dzi- 
bilchaltun is decreasing slightly (estimated value 
X — 0.91), and it seems that the seedlings are 
the most vulnerable stage of the life cycle, based 
on the observed mortality of this category. This 
information has to be considered carefully as de- 
mographic parameters (growth, survival, and fe- 
cundity) may change over time and space (Ca- 
swell 1989, Durán & Franco 1995). 
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